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Abstract 



A free fermion without doubler is formulated on 1+D dimensional discrete 

Minkowski space-time. The action is not hermitian but causes no harm. In 

1+3 dimensional massless case the equation describes a single species of Dirac 

particle in the continuous space-time limit. In 1+1 dimensional massless case 

the equation is the same as the automaton equation by 't Hooft and describes 

a chiral fermion. The time evolution operator is unitary and the norm is 

conserved. For interacting fermions with gauge fields the evolution operator 

is not unitary. If it is considered as an approximation for the theory on 

continuous space-time, the path integral formalism can be applied, where the 

fermion without doubling is used. Consequences of loosening the unitarity 

condition on the time evolution operator is discussed. 
11.15.Ha 
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I. INTRODUCTION 



As is well known, we necessarily meet the difficulty of so-called fermion doubling when we 
formulate a fermion on lattice space-time. Nielsen-Ninomiya proved under certain assump- 
tions that a chiral fermion cannot be formulated on Kogut-Susskind lattice (=continuous 
time and discrete space) At present it is widely believed that fermion doubling is in- 
evitable on discrete space-time in general. In order to get rid of this unwanted partner of 
the fermion Wilson introduced an additional term in Lagrangian to make the mass of the 
partner very large M. Hence the partner has no effect at low energy. Another attempt 
is one by Susskind f3j. He showed these two fermions could be interpreted as isodoublets 
(=fermions of two different flavours). 

In the present paper we want to show that a fermion can be formulated without doubling 
on Minkowski lattice space-time (Sec. 2). In order to avoid the no-go theorem we do not 
assume the hermiticity of action which was one of the assumptions of Nielsen-Ninomiya. 
The non-hermitian action may lead to some problems, but in our case it seems not so 
serious. For example, the action and the hermitian conjugate action yield two consistent 
field equations. This means both the real and imaginary parts of the action take stationary 
values at the same time, if the equation is satisfied. However, it is not yet clear whether the 
non-hermiticity of the action produces further difficulties. 

The equation we obtained on 1+D dimensional discrete Minkowski space-time has no 
fermion doubling solution. In 1+3 dimensions it tends to Dirac equation in the continuous 
limit. The time evolution operator is unitary and therefore, the norm is conserved. In 1+1 
dimensions a chiral fermion can be described in massless case which is the same as the 
automaton model of 't Hooft Q. However, in 1+3 dimensions a completely chiral fermion 
cannot be described as in the case of Susskind. This situation is not improved even if the 
unitarity of time evolution operator is not assumed (Sec. 4). 

In the above mentioned we treated only a free fermion field. In the case of interacting 
fields, however, the situation is completely different. It is difficult to assume the unitarity 
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of the time evolution operator. If we stand on the viewpoint that the discrete space-time is 
an approximation of continuous case, we may disregard this point, because the unitarity of 
time evolution operator is recovered in the continuous space-time limit. In the lattice gauge 
theory which is formulated on Euclidean space-time, the unitarity of time evolution operator 
is not considered seriously. Using the method of path integral formalism we can calculate 
any quantities in principle without the effect of fermion doubling, though the non-hermiticity 
might bring some practical difficulties (Sec. 3). 



In this section we develop a lattice theory of non-interacting fermions without doublers 
in 1 + D dimensional discrete space-time. We take a hyper-rectangular lattice with lattice 
constants r and a for the time and space direction, respectively. Lattice points are repre- 
sented by a set of 1 + D integers (tjX 1 , . . . ,x D ) = (t, x) = x. At each site is attached a 
spinor variable *&(x), whose number of components is not specified for the moment. 



where U is the time evolution operator acting on the x and spinor space. Postulating some 
properties on U, we will determine U and the minimal dimension of the spinor space for 

The time evolution operator U is assumed to be linear in the unit shift operators in the 
x-space so that only the nearest neighbor sites are coupled in the action, which will be given 
in the next subsection. Thus we have 



II. NON-INTERACTING FERMIONS 



A. Equations of motion 



We assume that equations of motion for \1/ take the form: 





D 




(2.2) 



i=i 



where the unit shift operators S^'s are defined by 
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(t, x\...,x\...,x D ) = tt(t, x\ . . . , x i + 1, . . . , x D ), (2.3) 

and Ai, Bi and C are matrices with respect to spinor indices. 

We furthermore assume that the time evolution operator is unitary: 

UlP = 1. (2.4) 

This condition is necessary for the theory to be well-defined quantum theory on discrete 
time. However, if the theory is considered to be an approximation to the continuum theory, 
we could loosen the condition in such a way that the unitarity should be recovered only in 
the continuum limit. In this section we assume the unitarity, and consequences of loosening 
the condition will be discussed in Sec. 4. 

We require that each component of \I/ satisfies the discrete version of the Klein-Gordon 
equation: 

1 /-2^- _£ S ,-2 + 4 +M ^ 0| (25) 

T 1=1 ° 

where M is the hermitian mass matrix. The dispersion relation implied by this equation is 
given by 

^ sin 2 ^ - -i f] sin 2 £ - M 2 = 0, (2.6) 
where k° and k are introduced through the Fourier transform: 

*(t,x)= ftfdke-'^ffcOk. (2.7) 

J — 7T 



Evidently there is no doubling problem with Eq. fl2.6|) , which has also been pointed out by 
Yamamoto ||. Therefore, it is this equation ( |2.5|) which assures that no doublers appear in 
our formalism. In the following we assume M = and at the end of this subsection we will 
discuss the case of M ^ 0. 

Substituting U given by Eq. (|2.2j) in the Klein-Gordon equation ( |2.5| ), where U~ l is 
replaced by U\ we find the relations: 
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At + B] = r, 
A + Bi=r, 

C + C ] - 2 = -2D, (2.8) 



where r is defined as (r/a) 2 . 

On the other hand the unitarity condition (|2.4j) requires that 

AiB] + AjB\ = 0, 

AiAl + BjBl = 0, (fori/j) 
AiCT 1 " + CB\ = 0, 

^(A4+ASt)+CCt = l. (2.9) 

i=l 

Combining those relations ( |2.8| ) and ( |2.9| ), we find that C/ takes the following form: 

U = E {^(1 + + £(1 - 4)4} + (1 - Dr) + iCj, (2.10) 

where the operators in the spinor space a«'s and Cj (the imaginary part of C) are yet to be 
determined by the following algebra: 

{ai,aj} = 2, 

{a l ,C I } = -2i(l-Dr), 

{oi,a]} = -2, (forz^j) 

f] {a,, 4} = D (- - 4D - 2) - ±Cl (2.11) 

i=1 \r / r 

Here a set of curly brackets { , } is an anticommutator. 

Using this algebra, we find that Cj is not linearly independent of a, but expressed as 

D 
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C/ = ^E( Q i-4)> (2-12) 
A 1=1 

which can be readily verified by showing the square of the difference between both sides 
of the equation vanishes. The relation fl2.12| ) has an interesting consequence on the lattice 



constants r and a. Calculating C] by Eq. (|2.12|) , one finds 



Cj = Dr(l- Dr), (2.13) 

which should be positive definite since Cj is hermitian. Thus we have 

r ^ < ±. (2.14) 

Yamamoto has obtained this inequality as well as the more general one mentioned below 
involving the finite mass ||. The argument is that the k° should be real in the dispersion 
relation ( |2.6[) , which is consistent with our requirement of the unitarity of U. 

It is convenient to decompose a, into the real and imaginary part as aj = Xi + iYi, where 
Xi and Yi are hermitian. By Eqs. fl2.12j ) and (|2.13|) the algebra for Xi and Yj is reduced to 
be 

{Xi, Xj} = 2- S^, 
{Xi,Yj} = 0, 

{y ij y i } = 2(i«s tf -i). (2.15) 

Although the anticommutation relation between X^s is diagonal with respect to the indices 
i and j, the one between Y^s is not. The matrix 2(5y/r — 1), however, can be easily 
diagonalized. The eigenvalues are 2(l/r — D) for the eigenvector (1, 1, ... , 1) and 2/r for 
any vectors orthogonal to (1, 1, . . . , 1). Thus the latter eigenvalue is (D — l)-fold degenerate. 

Therefore, when r < 1/D,we have 2D non-zero hermitian operators whose anticommu- 
tation relation is diagonal. Namely Xi and Yi, after the diagonalization, form the Clifford 
algebra up to trivial normalization factors: 



rj = r i , i = i,...,N 



t 

i 

{r i ,r i } = 2^ (2.i6) 

with iV being 2D in this case. So the dimension of the irreducible representation for X and 
Y is 2 D and *f>(x) should have 2 D components at least. 

In the case of r = 1/D, where r takes the maximal value allowed by Eq. ( 12.14] ), one 



of the eigenvalues 2(1 /r — D) vanishes and the linear relation = holds. In this case 

6 



we have the Clifford algebra with iV = 2D — 1, whose irreducible representation has the 
dimension of 2 D ~ 1 . Consequently ^>(x) has 2 D ~ 1 components. 

Assuming r — 1/D, so that *f?(x) has the smaller number of components 2 D ~ l , let us 
investigate the explicit form of the time evolution operator U for the cases D — 1 and D = 3. 

1+1 dimensional case 

x 1 ) is one-component field and U is given by 



where Xf = 1. When X\ = 1, we have U = Si and the equation of motion is simply 
+ 1, x l ) = ^>(t,x l + 1), which clearly simulates without doublers the chiral left moving 
fermion (8q— di)ty = in the 1+1 continuum theory. It is interesting that the correspondence 
is established by replacing the time derivative do^ by the forward difference quotient but 
the spatial derivative di^ by the backward one. This equation of motion has also been 
given by 't Hooft et al. in the cellular automaton approach [|J. In the same way, the case 
of Xi — — 1 describes the right moving chiral fermion without doublers. 

1+3 dimensional case 

The ratio of lattice constants r/cr = y/r is given by 1/Vo and ^(t, x) is 4-component 
field. The time evolution operator U then takes the form: 



U = ^(l + Xi)Si + ^(l-Xi)Sl 



(2.17) 



D 



Si s^ 



rVi 



(2.18) 



i=i 



i=l 



where l^'s are given by 




-Yi - Y 2 . 



(2.19) 



More explicitly those IYs (i 



1, . . . , 5) are expressed by the Dirac matrices: 
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r i = 7 V (« = 1,2,3), r 4 = 7 °, r 5 = 7 5 . (2.20) 

The second term in the time evolution operator fl2.18| ) is reduced to be Tj°j l di in the 
continuum limit. The last term in U goes to zero as r 2 in the naive continuum limit, and 
it reminds us of the Wilson term except for the spinor operators Yj in front. It is this 
spinor dependence that eliminates doublers completely, whereas the Wilson term only gives 
doublers a large mass. 

Now we briefly mention only the results for the case of M 7^ 0, which is more involved. 
Instead of the algebra fl2.ll]) , we have 



{ai,cij} = 2, 

{a i ,C I } = -2iC R , 

{oi,oJ} = -2, (forz^j) 

J2{a t ,4} = (l-C^C)- 2 -2D, (2.21) 

and the real part of C is given by Cr = (1 — Dr) — r 2 M 2 /2. An interesting consequence is 
that M 2 commutes with and Cj, and therefore, M 2 is just a number in the irreducible 
representation. We can furthermore show that r < (1 — r 2 M 2 /4)/D and even if the equality 
holds, ty(x) has 2 D components at least. 

B. Action 

We have shown that one can construct equations of motion for \l/(x) without doublers 



by the time evolution operator U of the form (2.2), which is local and unitary. Now we 
construct the action on the lattice which leads to the equations of motion. 

We immediately see that the desired equations of motion derive from the variation with 
respect to W of the action: 

S = iJ2 * f (t + k, x) (#(t + 1, x) - UV(t, x)) , (2.22) 



S 



where k is an arbitrary integer. The above action is evidently local but not hermitian. It 
should be noted, however, that the variation of the action with respect to \& also gives the 
equivalent equations of motion, which is a consequence of the unitarity of U. 

As we have seen in the last subsection, the single chiral fermion can safely reside in 1+1 
dimensional discrete space-time. This does not contradict the Nielsen-Ninomiya theorem, 
since our action is local but not hermitian. 

We will now show that the quantization through the path integral with the action (|2.22|) 



is equivalent to the equations of motion ( p.l|) with the canonical equal-time commutation re- 
lations {ty(t, x), x')} = 5 XX /. Here and in what follows, ^ and & should be understood 
to be canonically quantized fields whereas \l/ and W are Grassmannian variables. 

First we introduce the time evolution operator U in the fermion Fock space, which should 
not be confused with U, as follows: 

U = exp^-iJ2^(t^) H ^(t,^j , (2.23) 

where the hermitian operator H in the x and spinor space is defined through U = exp(— iH). 
The U is evidently unitary and reproduces the equations of motion: 

if(t + 1, x) = U-^t, x)U = UV(t, x). (2.24) 

We also find that U is independent of t. 

Now we calculate the matrix element of U between fermion coherent states | > and 
| >■ Here the coherent states are defined as 

| *(•) >= exp ^- Y, ^(x)^(0, x)^ | >, *(0,x)|0>=0. (2.25) 



The matrix element of U can be readily calculated as 



< *'(•) | U | > = exp (j2 ^' t (x)t/^(x) 



(2.26) 



Noting that the completeness relation of the coherent states is given by 

1 = /il [^(x)^(x)] exp U^^ t (x)^(x) J | *(•) >< |, (2.27) 



for the transition amplitude from the initial state | ^(to, •) > to the final state | ^(t n , ■) > 
we have 



exp (j2 x)tf(t n ,x) J | P^tp^r 



- 2 * f (* + !> x ) (*(* + 1, x) - C/*(t, x)) , (2.28) 

i=i ,x J 

which establishes the equivalence of the two ways of quantization after redefining W(t-\-k, x) 
as *t(t + i jX ) in the action (ggg) . 

Finally we will give the explicit form of the propagator, which is the inverse of the kernel 



in the action ( 2.22 ). Defining the unit shift operator in the time direction as So^(t, x) 
\l/(£ + 1, x), we have 

{ 1 \ 



< T$(x)¥(x') > 



i-slu, 



x,x' 

s -W 



(2.29) 



\(So-2 + Sl)-(U-2 + W)J XjX / 
Evidently the propagator has no extra poles, since in the momentum space its denominator 
takes the form of the left hand side of the dispersion relation (|2.6|), due to the Klein-Gordon 



equation (|2.5|). 



III. INTERACTION WITH GAUGE FIELDS 

We have seen in the previous section that the lattice fermion without doubling can be 
formulated for the free case in 1+D dimensions. In this section we consider the interaction 
of the fermion with gauge fields. 

The interaction is introduced by replacing the shift operators by covariant ones: 

S ^S (x)=V(x,x + 6)S , (3.1) 
S t (x) = V(x,x + t)S u (3.2) 
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where and % are unit vectors along the time and i-th space direction, respectively, and 
V(x, y) is a link variable connecting sites x and y. After these replacements the action is 
written as 

S = i J £tf( x ){l - Sl(x)U(x)}*(x) + S gange , (3.3) 

X 

where U(x) is defined by replacing Si with Si(x) in Eq. ( |2.2| ) and S^auge is the action for gauge 
fields. It is easily seen that the action is invariant under the lattice gauge transformation 
with a gauge group element g(x) : 

-> g(x)V(x), (3.4) 

V(x, y) -> (/(z) V(ar, y)g(y)~ 1 . (3.5) 

It is also checked that the new action has the correct continuum limit. 

Using this new action and path integral formalism, we can calculate the vacuum expec- 
tation value of an arbitrary operator O by the formula: 

<0>= J V^V^VV O exp(iS) I J V&VWVexp(iS). (3.6) 

In the usual case the gauge interaction is introduced in an analogous way using the 
symmetric discretization of time. But the unitarity of the time evolution of fermion fields is 
not clear for the symmetric discretization, since the field at time t + 1 depends on fields at t 
and t — 1. Usually this is not considered seriously because the unitarity is recovered in the 
continuum limit. However, in our case we have the explicit time evolution operator for a 
finite lattice spacing, and therefore, we can go a little further in the discussion of unitarity. 

We assume the Klein-Gordon condition and the unitarity condition also in this case. The 
Klein-Gordon condition gives the same relations Eq. (|2.8|) for A , B and C as in the free 



case. From the unitarity requirement we have the following relations : 

AiB]V(x, x + i)V(x + i,x + i + j) 

+A J B]V(x, x + j)V(x + j,x + i+j) = 0, (3.7) 
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AiAjV(x, x + i) V(x + i, x + i — j) 



+B j BjV(x, x - j)V(x -j, x + i-j) = 0, (i ^ j) 



(3.8) 



+ CB\ = 0, 



(3.9) 



D 



(3.10) 



i=l 



which correspond to Eq. Q2.9| ) in the free case. 



The last two relations are the same as before, but the first two contain link variables 
and depend on sites. The site dependence of the conditions means the unitarity of the time 
evolution operator U(x) cannot be satisfied by simple algebraic relations between A and B 
in general. 

1+1 dimensional case 

In 1+1 dimensions the relations reduce to the previous ones owing to the uni- 
dimensionality of space. The unitarity of the time evolution operator is kept by the same 
choice of A and B as in the free case. So the path integral formulation is well established 
keeping the unitarity at a finite lattice spacing. 

1+D dimensional case with D > 2 

We cannot determine A and B keeping the unitarity of the time evolution operator in 
this case. If we adopt the relations ( |2.9|) for A and B in the free case, we obtain from the 
two relations (3.7) and (3.8) 



[AiA] - BjBj){V(x, x + i)V(x + 1, x + % - j) 



V(x, x - j)V(x - j, x + i - j)} = 0, 



(i + j) 



(3.11) 
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(ABj - AjBj){V(x, x + i)V(x + i,x + i+j) 

-V(x, x + j)V(x +j,x + i+])} = 0. (3.12) 

We can see the origin of the unitarity violation comes from the difference between clock- 
wise and anti-clockwise parallel transports around a plaquette. This factor is higher order 
in the lattice spacing a and expected to be harmless in the continuum limit. 



IV. DISCUSSION ON UNITARITY 

In Sec. 2 we imposed the unitarity condition ( |2.4|) and the Klein-Gordon condition ( |2.5| 
on the time evolution operator U of the spinor field 

D 



U = Y,(A l S t + B l Sj) + C. (4.1) 
i=i 



We could not find the spinor field on discrete space-time with the same dimension as one 
on continuous space-time. 

In this section we try to find the spinor field with lower dimension than the one in Sec. 2 
without imposing the unitarity condition, since the Klein-Gordon condition ensures that the 
spinor field has no doubling. 

As in Sec. 2 we assume the Klein-Gordon equation: 

D 

U- 1 -2 + U = r Y,{Si ~ 2 + Sj) - r 2 M 2 , (4.2) 

i=l 

where we supposed the existence of inverse matrix for the time evolution operator. From 
the above equation we can see that the inverse matrix U~ l is of the same form as U : 

U-^'EiAiSi + Bity + C. (4.3) 
i=i 

Generally the term proportional to (S^) 71 (n > 1) may be included in £7 _1 , but these terms 
are not allowed, because in Eq. ( |4.2| ) there is no term which cancels out those terms. 
From Eq. (|4.2j ) we obtain the following conditions: 
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Ai + Ai 



Bt + Bi = r, 



C + C = K, 



(4.4) 



where K is 



K = 2(1 - rD) - t 2 M 2 



(4.5) 



and from UU 



-i 



1 we obtain 

AiAj + AjAi 
BiBj + BjBi 
~\~ C = 
BiC + CBi = 
AiBj + /^.l 



0. 



(for all i, j) 
(for all i, j) 
(for all i) 
(for all 2) 

(i + J) 



(4.6) 



We find that the matrix {C — K/2) 2 commutes with Ai, Bi and C, and we can regard 
this matrix as the unit matrix up to a constant factor DQ: 



(C-f) 2 = DQTL. 



(4.7) 



Using Eqs. 



> (|4.7| ) and assuming the isotropy of spatial lattice directions we can get 



the anticommutation relations between Ai, B^ and C 



{M a , M b } = F a 5 a>b 1, 



(4.8) 



where 



AD D 



(1 < a < D) 

2 (w-15-Q)i {D + l<a<2D-l) 
\+, (a = 2D) 



(4.9) 



(a = 2D + 1) 
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and A± are eigenvalues of the 2x2 matrix /A, which is defined by anticommutation relation 
between C - K/2 and E£i(A + Bj/y/D 



( 



M = 



/ 2 (L»r 2 -g+^) VDKr) 



V 



(4.10) 



'DKr 2L>Q J 

Clearly A± are real since the matrix M is hermitian . The matrices M a 's are defined by 

(1 < a < D) 



M n = I 



A a B a , 



Ef=iC-D(A + A), 



V? Ef =1 ? (A + A) " VDr + tf\C - f ), 



(2), 



V 1 ? Ef =1 r (^ + Si) - v^r + ^(C - f ), 



,( 2 )/ 



(D + 1 < a < 2D - 1) 
(a = 2D) 
(a = 2D + 1) 



where 77^ are the eigenvectors of M with eigenvalues A±: 



M 



( d)\ 



A 



± 



( <d\ 

1 T] K ± ' 1 



and the £^'s, (a = 0, 1, D — 1) are orthogonal vectors: 



We do not give the definite form of these vectors except for 



v 1 ; 



(4.11) 



(4.12) 



(4.13) 



(4.14) 
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since we do not need those forms in the later discussion. 
From these equations we obtain 

{MlMl} = F a 5 ab t. 



(4.15) 



There is no condition which prescribes the anticommutation relations between M a and M\ 
. For simplicity we assume the following anticommutation relation: 



{M a ,Ml} = D ab H. 



(4.16) 



Noticing that M a tends to a linear combination of 7 matrices in continuous space-time 
limit, we can easily see that the case with Q = (K 2 /4 — 1)/D is undesired. In this case all 
matrices M a vanish in the continuous space-time limit except for the two matrices M 2 d and 
M2D+1, although we need D + 1 matrices in order that the equation of spinor field goes to 
Dirac equation in the limit a, r — > 0. 

We consider the case with F a ^ (a — 1, 2D + 1). It is convenient to use 

M n 



M a = 



= X a + iY a . 

The anticommutation relations between these matrices become 



^ ^ab + Gab Fab ^ 





[.0 


} 


[ 


,(x b ,Y b ) 









ab Vab 



(4.17) 



(4.18) 



where 



and D ab is 



Gab = -^ab ~ D ab - D* ba )t, 
F a b = ^(D ab -D* ba )ti, 



D = Dab 

a VKF b 



(4.19) 



(4.20) 
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As the eigenvalues of matrix Gab are non-negative, the matrix 8 a b + Gab is invertible. Using 
the inverse matrix of 5 a b + Gab, we define 



6=1 



(4.21) 



and then we have 









( 


[ 




, {x t , n) J = 






{ K ) 




V 



/ 



(4.22) 



5ab + Gab 

Vab + (^(1 + GY X T)ab 

As the upper-left submatrix of the right hand side of the above equation is invertible, we 
cannot make the number of independent matrices X a and Y a smaller than 2D+1. Thus the 
dimension of these matrices is larger than 2 D ( for example, the lower limit of the dimension 
is 8 for D = 3). 

If A + and A_ are zero, the dimension of matrices is equal to or larger than 2 D ~ 1 . If either 
A + or A_ is zero, the dimension of matrices is equal to or larger than 2 D . For example, the 
limit of the dimension is 4 or 8 for D = 3. The massless spinor field in Sec. 2 corresponds to 
the former case. 

In conclusion we found under the assumption ( j4.16| ) that the dimension of spinor field 



without the unitarity condition is not smaller than the dimension with the condition. 



V. SUMMARY 

We have formulated a free fermion without doubling on 1+D dimensional Minkowski 
lattice space-time. We required there the unitarity of time evolution operator and Klein- 
Gordon equation on lattice space-time. This means the norm is conserved and the fermion 
has no doubler. We showed that the minimal number of components of massless field ^ is 
2 D ^ 1 . In 1+1 dimensional case the equation is the same as that of the cellular automaton by 
't Hooft et al. In 1+3 dimensional massless case the time evolution operator was expressed 
in an explicit form using usual 7 matrices, which tends to Dirac operator in the continuous 
space-time limit. 
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The action is not hermitian. We proved the equivalence of canonical quantization and 
that through the path integral. We have given the explicit form of the fermion propagator, 
which has no extra poles of doublers. 

In the case where the fermion interacts with gauge fields the action was also written in a 
gauge invariant form. The vacuum expectation value of an arbitrary operator is calculated 
by the path integral formulation in principle. The time evolution operator is not unitary in 
this case, and the reason was considered briefly. 

We have tried to find the spinor field with lower dimension than without imposing 
the unitarity condition. However, we found the spinor field with the smallest dimension is 
2 D ~ 1 under a certain condition. 
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